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Abstract 

We introduce a probabilistic representation for solutions of quasi- 
linear wave equation with analytic nonlinearities. We use stochastic 
cascades to prove existence and uniqueness of the solution. 
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We consider the following nonlinear wave equation on the real line R: 
Uu(x, t) = F(x, t, u(x, t)), (x, t) £Rx R + . (1) 

Here 

Uu(x, t) = u tt (x, t) - u xx (x, t), 

and F is a given function. 

For T > 0, we say that u : M. x [0, T) is a classical solution of the Cauchy 
problem of ([I]) with initial conditions 

u(x,0) = 4>{x), iGi, (2) 
u t (x,0) = tp{x), xER, (3) 

if u G C 2 (M x (0, T)) n C l (R x (0, T)). 

The goal of this note is to develop a stochastic cascade approach to 
constructing solutions of the Cauchy problem dl])-©. It is similar to the 
construction of solutions for the Navier-Stokes system suggested in [LJS97] 
and for the KPP equation in [M75j . Although this approach is essentially 
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equivalent to a Picard-type iteration scheme, it provides an interesting point 
of view. 

Probabilists have long desired a probabilistic representation of the wave 
equation, but there are only a few papers on the topic. In |Kac56j IKac74] , 
Kac discovered a probabilistic representation for the telegrapher's equation, 
which is a wave equation with a lower-order time derivative. More recently, 
Dalang, Tribe, and the second author developed a multiplicative version 
of the Feynman-Kac formula which applies to the wave equation, among 
others, see [DM091 IDMT08] . 

We begin with the classical d'Alembert representation of solutions for 
the linear wave equation. If F = 0, i.e., the problem ([l])-(j3]) is a homoge- 
neous linear problem, and existence and uniqueness hold under unrestrictive 
assumptions on the regularity of initial conditions (j) and tp. Fixing (j> and 
■0, one can write the solution of the homogeneous Cauchy problem as 

v{x,t) = \! il){y)dy + \{ct>{x + t) + (l>{x-t)), (4) 

1 Jx-t 1 

If F(x,t,u) = f(x,t) is a sufficiently smooth function that does not 
depend on u, then we have an inhomogeneous wave equation with external 
source /, and the d'Alembert formula holds: 

u(x,t) = v(x,t) -\ — / f(x,t)dxdt, (5) 

2 JA(x,t) 

where 

A(x, t) = {(y, s) : < s < t, \y -x\ <t- s} 

is the light cone of the past associated with the space-time point (x,t). In 
fact, in our one-dimensional situation, A(x,t) is just a triangle. 

Formula ([5]) allows us to define a mild solution of equation jl} ona time 
interval [0, T) as a measurable function u : M. x [0, T) — > R such that for all 
(x,t) E R x [0,T), 

u(x,t) = v(x,t) + - f F(x,t,u(x,t))dxdt. (6) 
2 J A(x,t) 

From now on we shall assume for simplicity that F(x,t,u) = F(u) does 
not depend on (x,t), although our construction can be also applied with 
appropriate modifications in the general case. The next assumption is crucial 
for our construction, though: we require analyticity of F, i.e., we assume 
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that for all it, F(u) can be represented as a convergent power series in u: 

oo 
k=0 

In particular, we can deal with power-type nonlinearities like F(x, t, u) = u 2 . 

Let us fix (x, t) £ lxR + and try to rewrite the d'Alembert formula in the 
language of random variables. To that end, let us introduce a random point 
(£, r) uniformly distributed in A(x, t). Since the area (Lebesgue measure) of 
A(x, t) equals t 2 , it means that the random point (£, r) has density vanishing 
outside of A(x,t) and identically equal to t~ 2 inside A(t,x). Therefore, 
formula ([6]) can be rewritten as 

t 2 

u(x,t) = v(x,t) + -EF(u(C,T)). (7) 

Next step is to consider a sequence of numbers (pk)k*Lo with the following 
properties: 

(i) p is a probability distribution: Ylk^oPk = 1> Pk > for all > 0; 

(ii) pfc > for every k > with a/c 7^ 0; 

(iii) J2T=i k Pk < 1- 
Also let us define 

w ( x t) = v( - x ' t \ f x t) e R x R+. 

and, for each k > 0, 



0, = 0, 



Let us introduce a random variable k distributed according to p and 
independent of (£, r). Then ([7]) immediately implies the following lemma. 



Lemma 1 If u is a solution of ((TJ) on [0, T), then for any (x,t) 6 



it(x, i) 



(8) 
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Next step in our construction is to iterate Lemma [TJ Namely, for any 
k > 0, on the event {k = k} we may compute r) by the same procedure. 
The role of (x, t) is played by (£, r), and, conditioned on {(£, r) = [x 1 ', t')}, to 
compute the product of k copies of u(x',t') we may consider k independent 
random variables Tj, /c$)* =1 , so that random points (&,7i) are uniformly 
distributed in A x / ^ , and random variables distributed according to 

distribution p. Given that collection of random variables, for each i = 
l,...,k, we can apply Lemma [TJ Notice that on the event {k = 0}, the 
random variable under the expectation sign in ([8]) is a constant equaling 
w(x, t) + L ^-, so that we do not have to consider any new random variables 
to compute it. 

It is clear that iterating this procedure we obtain a stochastic cascade, 
i.e., a branching process with each particle assigned a location in space and 
time. To make this idea precise, let us introduce more notation. We shall 
need a probability space rich enough to support these random structures 
involving random family trees of the participating particles and their random 
locations. 

We begin with an encoding of vertices of finite rooted ordered trees. 
Each vertex v in the n-th generation of the tree can be identified with a 
sequence (a±, . . . , a n ), where cij G {0, 1, . . .} for all i = 1, . . . , n. The parent 
of (ai, . . . , a n ) is (a%, . . . , a n -\). It is convenient to identify the root of the 
tree with an empty sequence denoted by which is consistent with the above 
encoding of the parent-child relation. 

Let us fix n £ N and let probability measure on rooted ordered trees with 
at most n generations be given by the classical Galton- Watson distribution 
on trees with branching distribution (pfc)j£L > i- e -> eacn vertex v of the tree 
has a random number k v of children (v, 1), . . . , (v, k), where (v, i) means a 
sequence obtained from v by appending (concatenating) i on the right. The 
random variables k v are jointly independent. 

Each vertex v in thus generated random tree gets a random space- 
time label (k v ,t v ) according to the following rule. First, we set (kq,tq) = 
(x,t). Then, we can iteratively apply the following: for any vertex v = 
(ai, . . . , a m ), m < n, conditioned on k v ,t v = (y, s) and on k v = k, the labels 
T (v,i))i=i are i.i-d. uniform random points in A(y, s), independent of 
all previously constructed random elements in the procedure. 

Now we shall recursively define a functional IT on subtrees of individual 
vertices starting with the leaves of the tree. There are two types of the 
leaves. Leaves of type 1 are vertices of generation n. For any leaf v of 
type 1, we set H(v) = u(£ v ,t v ). A leave of type 2 is a vertex v that did not 
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produce any children, = 0. For these vertices we set 

U(v)=w^ v ,r v ) + ^. (9) 

As soon as the values of H(v) have been assigned for all leaves of the realiza- 
tion of the tree, we may start assigning values to other vertices iteratively. 
For any vertex v of the tree such that H(v, i) has already been assigned for 
all i = 1, . . . , k v , we define 

nW = $ K „n«K H ,r M ). (io) 
z i=i 

Proceeding iteratively, we eventually will assign some value 11(0) to the root 
of the tree. This value is a random variable incorporating information from 
all the other vertices of the random tree as well as their space-time labels. 
We denote the resulting random variable by H n (x,t,u(-)) to stress that the 
tree of at maximum n generations was initiated at (x, t) and the solution u 
was used to evaluate II at the tree's leaves. 

Theorem 1 Suppose u is a solution of (pQ) on [0, T] . Then for any n and 
any (x,t) £lx (0,T), 

u(x,t) = EU n (x,t,u(-)). 

An obvious next step is to take n to infinity. Since the requirement (iii) 
on the branching distribution p means that the branching process is critical 
or subcritical, the realizations of the random trees almost surely have finitely 
many vertices. In particular, with probability 1, 

lim U n (x,t,u(-)) = U OQ (x,t), (11) 
n—*oo 

where the random variable H OQ (x,t) is constructed from the realization of 
the stochastic cascade in exactly the same way as H n (x,t,u(-)) for finite n 
except that there are no leaves of type 1. 
If we show that 

lim EII n (x,t,u(-)) = EIIooOM), (12) 

n— too 

then we will be able to conclude that 

u(x,t) = En oo (x,t). (13) 



5 



In that case, since II^Oe, t) is a (random) functional of w, the (modified) 
external source only and does not involve u, we can claim that the solution 
is unique and it is given by formula (|13p . It is also easy to see from ([8]) that 
if (fl2|) holds then u(x, t) given by formula (fl~3j) is a solution. 

So, condition (fT2"j) implies the existence and uniqueness of solution and 
its stochastic representation. 

We notice that P{n oo ^ Il n } — > as n — > oo. Therefore, to ensure 
convergence in (|12p it is sufficient to check that £1100(2;, t) is well-defined. 

Next, we notice that, according to ([9]) and (fTOj) Ii{x,t) is a product of 
many factors of the form w(£ v ,t v ) + or T « . If we require that all 
these factors are bounded by 1, then n oo (x, t) product is bounded by 1, and 
its expectation is also bounded. This leads us to the following theorem. 

Let 6* = sup fcgN Let 

T* = Bupjt<y|: supKx,t)| + ^<l|. 

Theorem 2 IfT*> 0, then there is a unique solution of JT]) on [0, T*). It 
is given by (fT3|) . 

The crude requirement of boundedness by 1 and the resulting condition 
of the above theorem can certainly be improved for some specific cases. 
We do not explore this issue further since we do not expect our method to 
produce sharp conditions for the existence of the solution. 
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